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Mild traumatic brain injury (mTBI) has become an increasing public health concern as subsequent injuries can 
exacerbate existing neuropathology and result in neurological deficits. This study investigated the temporal 
development of cortical lesions using magnetic resonance imaging (MRI) to assess two mTBIs delivered to op- 
posite cortical hemispheres. The controlled cortical impact model was used to produce an initial mTBI on the 
right cortex followed by a second injury induced on the left cortex at 3 (rmTBI 3d) or 7 (rmTBI 7d) days later. 
Histogram analysis was combined with a novel semi-automated computational approach to perform a voxel-wise 
examination of extravascular blood and edema volumes within the lesion. Examination of lesion volume Id post 
last injury revealed increased tissue abnormalities within rmTBI 7d animals compared to other groups, particularly 
at the site of the second impact. Histogram analysis of lesion T2 values suggested increased edematous tissue within 
the rmTBI 3d group and elevated blood deposition in the rm TBI 7d animals. Further quantification of lesion compo- 
sition for blood and edema containing voxels supported our histogram findings, with increased edema at the site of 
second impact in rmTBI 3d animals and elevated blood deposition in the rmTBI 7d group at the site of the first injury. 
Histological measurements revealed spatial overlap of regions containing blood deposition and microglial activation 
within the cortices of all animals. In conclusion, our findings suggest that there is a window of tissue vulnerability 
where a second distant mTBI, induced 7d after an initial injury, exacerbates tissue abnormalities consistent with 
hemorrhagic progression. 

© 2012 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

Traumatic brain injury (TBI) affects an estimated 1.7 million indi- 
viduals and contributes to a third of injury-related deaths annually in 
the United States, with mild TBI (mTBI) accounting for approximately 
75% of all TBI cases (Corrigan et al., 2010; Faul et al., 2010; Tagliaferri 
et al., 2006). Comparison of clinical computed tomography (CT) and 
magnetic resonance imaging (MRI) scans found that MRI allows for 
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increased detection of mTBI related abnormalities (Lee et al., 2008). 
However given the standardized definition of mTBI wherein there 
are no observable changes on standard MRI, new approaches are 
warranted to assess mTBI from neuroimaging data (VA/DoD, 2009). 

Standard MRI, such as Tl and T2 weighted imaging (T1WI, T2W1), of 
mTBI results in low contrast subtle tissue changes that create difficulty 
in evaluation and detection of mild injury. As a result, enhanced MRI 
quantification methods are needed to identify and monitor the progres- 
sion of mild injuries. Currently, quantitative semi-automatic MRI analy- 
sis of moderate to severe experimental TBI has been used to investigate 
tissue pathology (Immonen et al., 2009; Irimia et al., 201 1 ; Kharatishvili 
et al„ 2009). In these experimental models, manually defined regions of 
interest revealed positive correlations between T2 relaxation values and 
poor histological and behavioral outcomes (Immonen et al., 2009; 
Kharatishvili et al., 2009). However, little work has been done using 
models of mTBI. One study in which an automated method of whole 
brain MRI analysis was developed, demonstrated increased abnormal 
gray matter within the first 5 h following experimental mild and severe 
TBI (Colgan et al., 2010). However, long-term (days, weeks) quantitative 
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computational MRI analysis of the evolution of neuropathology follow- 
ing mild injuries, and in particular repeated injury, remains unexplored. 

One of the many consequences of TBI is microvascular dysfunction, 
which can result in tissue ischemia, edema and hemorrhagic progression 
(Kurland et al., 2012). Depending on the type and severity of brain injury, 
peak edema is typically observed 24-48 h following injury and resolves 
by approximately 4-7 days, while hemorrhagic lesions can occur within 
the first 3-4 days and resolution can take months (Kurland et al., 2012; 
Obenaus et al., 2007; Oehmichen et al., 2003). Hemorrhagic progression 
of a contusion (HPC) is the expansion of an existing or the development 
of a new hemorrhagic lesion (Simard et al., 2010). HPC is observed clin- 
ically in mild to severe cases of TBI and substantially increases tissue 
damage, as blood is particularly neurotoxic (Alahmadi et al., 2010; 
Kurland et al., 2012; Sifri et al., 2004, 2006; Simard et al., 2010; Xue 
and Del Bigio, 2000). Experimental modeling has demonstrated that ves- 
sel fragmentation plays a key role in HPC development and blocking 
these cellular cascades results in decreased necrotic tissue volume and 
improved behavioral outcomes (Simard et al., 2009). 

Most patients who sustain an mTBI appear to recover normally 
from the injury, though there is evidence that sustaining multiple 
mTBIs can result in cumulative detrimental effects (Belanger et al., 
2010). However little is known about the underlying pathology and 
long-term consequences of repeated mTBI (rmTBI). Animal models 
where induction of two injuries to the same location was performed 
demonstrated a window of vulnerability of approximately 1-5 days 
following an mTBI, in which a second impact can result in increased 
lesion volume, axonal injury and inflammation (Laurer et al., 2001; 
Longhi et al., 2005; Shitaka et al., 2011). However two mTBIs 7 days 
apart did not result in exacerbation of injury (Longhi et al., 2005). Ad- 
ditional models of rmTBI in which four or five injuries to the same lo- 
cation were rapidly induced (24 h apart) demonstrated persistent 
neurological deficits with no overt cell death or blood brain barrier 
disruption (DeFord et al., 2002; Kane et al., 2012; Shultz et al., 
2012). Surprisingly, all of the reported models of rmTBI had repeated 
injuries to the same anatomical location and what remains to be elu- 
cidated are the effects of multiple injuries to different brain locations. 

The objective of our study was to evaluate whether the ongoing neu- 
ropathology from an initial mTBI was altered by a second contralateral 
injury induced 3 or 7 days later. A secondary objective was to assess 
the ability of computational analysis of the mTBI injury to provide 
quantitative information about the progression of neuropathology. We 
hypothesized that injuries induced 3d apart would result in larger lesion 
volumes at the site of the second injury, while little or no change in le- 
sion volume was expected within groups receiving impacts 7d apart. 

2. Methods 

Protocols were in compliance with federal regulations and approved 
by the Loma Linda University Animal Health and Safety Committees. 

2.3. Mild traumatic brain injury 

Fifty adult male Sprague Dawley rats ages 2-4 months were random- 
ly assigned to three experimental groups: Single mTBI and repeated 
mTBI induced either 3 (rmTBI 3d) or 7 (rmTBI 7d) days apart (Fig. 1A). 
A mild controlled cortical impact (mCCI) was used to induce mTBI as 
previously described (Obenaus et al., 2007), with minor modifications. 
Briefly, anesthetized rats (Isoflurane: 3% induction, 1-2% maintenance) 
were secured into a stereotactic frame and a midline incision exposed 
the skull surface. A craniectomy (5 mm diameter) was performed over 
the right hemisphere 3 mm posterior and 3 mm lateral from Bregma 
(Fig. IB), where a mCCI (4 mm diameter tip, 0.5 mm depth, 6.0 m/s 
speed, 200 ms dwell) was delivered to the cortical surface using an elec- 
tromagnetically driven piston (Leica Biosystems Inc., Richmond, IL). An- 
imals within the rmTBI groups received a second craniectomy and mCCI 
on the opposite (left) cortical surface (identical Bregma coordinates) 3 or 



7 days after the first mCCI (Fig. 1 A, B). At the end of each surgery a sterile 
nylon flap (-1 cm 2 ) was placed over the craniectomy and the incision 
was closed. 

2.2. MRI data collection 

Animals were anesthetized (isoflurane: 3% induction, 1% mainte- 
nance) for MRI at: Id post first injury (Dl ), Id post last injury (Dl Singles, 
D4 rmTBI 3d, D8 rmTBI 7d), 14d post injuiy (D14: 14d post first, D21 : 14d 
post last injuiy) and >60 (60-68) days post first injuiy (Fig. 1A). Due to 
imaging constraints, all MRI was performed within 24 h of the targeted 
time point with the exception of the rmTBI 3d group, which received im- 
aging at 17d rather than 21d post initial injuiy. In vivo T2 weighted im- 
ages (T2WI; TR/TE/FA = 3453 ms/20 ms/20°, 25x1 mm slices) and 
susceptibility weighted images (SWI; TR/TE/FA = 39 ms/20 ms/20°, 
48x0.8 mm slices) were acquired with a 256x256 matrix using a 
4.7 T Bruker Avance. SWI was used to identify extravascular blood within 
tissues and T2WIs were used to create T2 maps computed using in-house 
software (Obenaus et al., 2007). 

2.3. MRI data analysis 

Regions of interest (ROI), manually drawn on T2WIs using Cheshire 
image processing software (Hayden Image/Processing Group, Waltham, 
MA), included the right and left hemispheres and cortical areas 
containing observable abnormalities (hyper/hypo-intense signals). In 
a subset of animals (n = 20) transient edema was observed within the 
corpus callosum following the initial impact that did not expand follow- 
ing a second impact and quickly resolved; in this group we only 
assessed cortical abnormalities. For voxel-wise computational analysis 
of the lesions, both the ROIs and T2 maps were imported into a novel 
computational algorithm coded in MATLAB (MathWorks, Torrance, 
CA) where lesion T2 values, tissue level histograms and volumes were 
extracted and lesion composition was classified (Fig. 1C). 

T2 value ranges for blood, edema and normal appearing brain ( NAB ) 
containing voxels were determined by taking the average T2 value 
of manually segregated blood and edema lesion pixels. This was 
performed using a training-set composed of MRI scans (6 MRI scans 
evaluated) collected at Id post last injuiy from animals in each group 
that contained lesion volumes representative of the group average 
and contained both edema and blood. From this analysis, lesion voxels 
were classified as blood (<69 ms), edema (>90 ms), NAB (70-89 ms) 
or noise ( > 500 ms). Following voxel classification, a connected compo- 
nent analysis was performed on each slice of the brain volume in which 
neighboring voxels of the same class (blood, edema, NAB, and noise) 
were grouped (Fig. 1C). To exclude voxels with a high spatial frequency, 
groups were then filtered by size such that any group less than two 
voxels were reclassified as noise (Fig. 1C,D). 

T2 histograms of the entire lesion volume were created (0-500 ms, 
binned into 1 ms bins) and kurtosis, skewness, mean and peak values 
were also calculated for each injuiy histogram. Total (first + second) his- 
tograms were analyzed at Id post last injuiy and 14d post first injuiy; 
first and second injuiy histograms were analyzed at Id post last injuiy 
and 14d post injuiy (first: 14d post first, second: 14d post last injuiy). 
Group histograms were normalized to the average number of pixels con- 
tributing to the lesion in order to account for differences in lesion size. A 
threshold of 150 pixels (approximately 2 ul of tissue) was used to reduce 
unnecessaiy noise within the data set by including only voxels clearly 
containing injured tissue. Histograms then underwent 2D smoothing 
using a low pass Gaussian filter to remove high frequency noise. 

Lesion, edema and blood volumes were divided by the whole 
brain volume to account for differences in brain size between animals. 
Total lesion volumes (first + second) were analyzed from animals at 
Id post last injuiy, 14d post first injury and at the time of final MRI 
(Single, rmTBI 7d: 60-68d and rmTBI 3d: 21 d post first injury). First 
and second lesion, blood and edema volumes were analyzed at Id 
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Fig. 1. Experimental design and schematic of computational analysis algorithm. A. Experimental mTBI and neuroimaging timeline. An mTBl was induced to the right cortex on DO 
(denoted as an *) in all animals. The second mTBI was induced to the left cortex at either 3 or 7 days later (*). MR1 was performed Id post first (Dl ), Id post last (Dl Singles, D4 
rmTBI 3d. D8 rmTBI 7d) and 14d post injury (D14: 14d post first, D21 : 14d post last injury) (red circles). B. Illustration of the mTBI locations for the first (right) and second (left) 
injuries. C. Computational schema forT2 MR1 analysis to assess lesion volumes and composition for blood (<69 ms), edema (>90 ms), normal appearing brain (NAB; 70-89 ms), or 
noise (>500 ms) containing voxels. D. RawT2 weighted MRI after a Single mTBI revealed observable abnormal signal intensities (yellow arrows) that were then subjected to com- 
putational analysis for tissue classification (edema, blood, NAB, noise). 



post last injury and 14d post injury (first: 14d post first, second: 14d 
post last injury). Calculations for total (first + second), first and second 
lesion volumes and blood and edema volumes were performed where: 
a) total lesion volume was calculated by dividing the total lesion volume 
(first + second) by the total brain volume, b) first and second lesion vol- 
umes were calculated by dividing the first or second lesion volume by 
total brain volume and, c) blood and edema volumes were calculated 
by dividing the volume of blood or edema found in either the first or 
second lesion and dividing this by the respective first or second lesion 
volume. 3D reconstructions of lesion volumes within the brain were 
created using Amira software (Mercury Computer Systems, Visage 
Imaging, Inc., San Diego, CA, USA). 

Animals with a total lesion volume (48.1 -58.5 ul) greater than 4% of 
whole brain volume at any time point were excluded ( n = 5 ) as this le- 
sion size fell two standard deviations away from the average total lesion 
volume Id following the initial injury and was therefore not considered 
a mild injury. Additionally one animal within the Single group demon- 
strated an increasing lesion volume over time, which was not seen in 
any other animals within the group and was also excluded from the 
study. While our study compared mild Single and repeated TBI animals, 



a cohort of Sham animals ( n = 3 ) was assessed undergoing the identical 
experimental procedure as Single animals but without TBI. Using our 
computational analysis on Sham cortical ROIs, we found no significant 
differences (p>0.05; p = 0.82-0.95) in blood, edema, NAB or noise vol- 
umes (ul) between the left and right cortical regions. 

2.4. Tissue collection 

Following the final MRI session animals were sacrificed via 
transcardial perfusion. Extracted brains were prepared for staining 
as previously described (Obenaus et al., 2007). Briefly, brains were 
post-fixed in 4% PFA overnight, undergoing two 30 min PBS washes 
following fixation. Prior to cryosectioning, brains were placed in 30% 
sucrose and embedded in optimal cutting temperature compound 
(O.T.C., Tissue Tek; Sakura Fine Tek, Torrance, CA). Coronal sections 
(30 nm) were collected using a Leica CM1850 cryostat (Leica 
Microsystems GmbH, Wetzlar Germany) and mounted directly on 
gelatin-chrome-alum-coated slides for histological staining (stored 
at —20 °C) or were free floating (stored in cryoprotectant solution 
at 4 °C) for immunohistochemistry. 



V. Donovan et al. I Neurolmage: Clinical 1 (2012) 18-28 



21 



2.5. Histochemistry 

All staining was performed on slices centered under the lesion site at 
14d post first injury. Prussian blue staining was performed as previously 
described (Obenaus et al., 2011) to detect extravascular blood within 
the tissue. Degenerating neurons were detected using Fluoro-Jade B 
(Millipore, Temecula, CA) staining on thawed tissue sections as previ- 
ously described (Schmued and Hopkins, 2000). 

Immunohistochemistry was performed on free floating sections as 
previously described (Obenaus et al., 2007). Immunolabeling was 
performed for astrocytes, using glial fibrillary acidic protein (GFAP; 
Millipore, Temecula, CA, 1 :1 000), and microglia, using ionized calcium 
binding adaptor molecule 1 (IBA1 ; Wako, Carpinteria, CA, 1 :400). Sec- 
ondary antibodies used included goat anti-mouse Alexafluor 488 
(Invitrogen, Carlsbad, CA, 1:400) and goat anti-rabbit rhodamine 
(Millipore, Temecula, CA, 1:200). Negative controls omitted the pri- 
mary antibody during the staining procedure. 

The extent of blood deposition (Prussian blue) and activated mi- 
croglia (IBA1 ) within the tissue was quantified using N1S Elements Soft- 
ware (Nikon Instruments Inc., Melville, NY, USA). A total of three 
animals (2 sections/animal), containing lesion volumes representative 
of the group average, were selected from each experimental group 
and were used to measure the medial-lateral and dorsal-ventral dis- 
persal of blood deposition and activated microglia within the injured 
cortices. Activated microglia were identified based on visual signal in- 
tensity and cellular morphology as previously described (Obenaus 
et al., 2008). 

2.6. Statistics 

To correct for variability following the initial mTBl, data points were 
normalized to each animal's initial recorded value (current value/initial 
value), such that all animals received a score of 1 for the first impact. A 
non-parametric Kruskal-Wallis analysis of variance (ANOVA) was used 
for statistical analysis of the total (first + second) lesion volume as this 
data failed to meet the requirements of normality using a Shapiro- Wilk 
test. A one-tailed t-test was used to test the rmTBl 7d second lesion 



volume against the Single first lesion volume, as we hypothesized that 
the site of the second impact would have a greater injury volume than 
a Single injury alone. However, two-tailed t-tests were used to compare 
the first and second lesion volumes of all groups at 14d post injury. Fur- 
ther, two-tailed t-tests were used to evaluate the first and second blood 
and edema volumes at Id post last injuiy. A p value < 0.05 was consid- 
ered significant (n = 22 Single, 4 rmTBI 3d, and 15 rmTBI 7d); error 
bars represent the mean ± SEM. 

3. Results 

3.1. Lesion composition is dependent on the interval between rmTBI 

Analysis of the mTBl lesion was first evaluated using total lesion 
(first + second) T2 value histograms to ascertain temporal changes 
and identify rmTBI group differences at Id post last injury. Compared 
to Singles, the rmTBI 3d animals exhibited a rightward shift in the T2 
total lesion histograms suggesting increased edematous tissue that 
was not seen in the rmTBI 7d group (Fig. 2, top panel). In contrast, a left- 
ward shift was observed in the rmTBI 7d lesion compared to Single an- 
imals (Fig. 2, top panel) consistent with increased blood deposition. We 
then further assessed each injury site independently and found that the 
T2 histogram at the site of the initial injuiy of both rmTBI groups had a 
leftward shift (increased blood) compared to animals receiving only a 
Single mTBl. However at the site of second injury, the rmTBI 3d group 
exhibited a dramatic shift rightward indicating increased edema, 
while the rmTBI 7d animals demonstrated a similar T2 profile to that 
observed in Single animals (Fig. 2, top panel). 

At 14d post first mTBl, total brain histograms revealed a strong 
leftward shift in all experimental groups, demonstrating that at this 
sub-acute time point the edema seen at Id post last injuiy had re- 
solved. Further, the strong leftward shift of the histogram (T2 peak 
values: Single = 52.6 ms, rmTBI 3d = 47.8 ms, rmTBI 7d = 55.5 ms; 
see Table 1 ) in all groups suggests that extravascular blood deposition 
substantially contributed to the overall lesion composition at this 
later time point (Fig. 2, bottom panel). A closer examination of the 
lesion created by the initial impact revealed that all groups exhibited 
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Fig. 2. Temporal shifts in T2 value distribution following rmTBI. Histograms of the total lesion (first+ second) at Id post last and 14d post injury (14d post first injuiy) demonstrate 
global changes in T2 value distribution between groups. The total lesion was then separated into the first and second lesions at Id post last injury and 14d post injury (first: 14d post 
first, second:14d post last injury) to illustrate the differences in T2 distribution between the two injuiy sites. Color coding in the total lesion histograms illustrates the range of T2 
values for tissue classification: blood (green), edema (red), or normal appearing brain (NAB; blue). Single mTBl histograms in the second injuiy graphs (dotted line) are those taken 
from the injured cortex (first injury) for comparison. 
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Table 1 

Total lesion histogram features. Peak T2 values, mean T2 values, skewness and kurtosis 
measurements collected from total (first + second) lesion within Single, rmTBI 3d and 
rmTBI 7d groups Id post last and 14d post first injury. Measurements ± SEM. 





Time point 


Single 


rmTBI 3d 


rmTBI 7d 


p Value 


Peak T2 value 


ID post last 


70.5 ±1.5 


74.2 ±6.1 


68.3 ±2.5 


p = 0.505 


(ms) 


injury 

14d post first 


52.6 ±3.0 


47.8 ±5.4 


55.5 ±2.8 


p = 0.342 


Mean T2 value 


injury 

ID post last 


85.0 ±1.4 


92.6±5.3 a 


78.8±2.6 a 


p = 0.015 


(ms) 


injury 

14d post first 


63.9 ±2.9 


55.8 ±2.9 


62.6 ±3.0 


p = 0.263 


Skewness 


injury 

ID post last 


2.3±0.1 a 


2.3 ±0.3 


2.9±0.2 a 


p = 0.042 




injury 

14d post first 


2.4±0.2 a 


5.3±0.6 a 


3.8 ±0.7 


p = 0.026 


Kurtosis 


injury 

ID post last 


14.8±1.2 a 


13.7 ±2.2 


21.3±2.6 a 


p = 0.033 




injury 

14d post first 


14.7 ±1.8"' 


47.8±10.4 a 


29.6 ±7.6 


p = 0.034 




injury 











a Designates measurements that are statistically different from each other. 



a similar T2 profile at 14d post first injury (Fig. 2, bottom panel). However 
at the second impact site, the rmTBI 3d animals demonstrated a leftward 
shift compared to its Id post last injury data and the rmTBI 7d group at 
14d post second injury, indicating the presence of blood and resolution 
of acute edema (Fig. 2, bottom panel). The site of the second impact in 
rmTBI 7d animals shared a similar T2 distribution to that of Singles 
(Fig. 2, bottom panel). 

Quantitative evaluation of total lesion histogram features was un- 
dertaken (Table 1 ) with no significant differences in the T2 peak values 
between any group or time point. However, the mean histogram T2 
values within the total lesion at Id post last injury revealed a significant 
(p = 0.01 5) increase in T2 values within the rmTBI 3d group compared 
to the rmTBI 7d group (Table 1 ). At 14d post first injury there were no 
significant differences between mTBI groups. Analysis of skewness 
( measure of distribution asymmetry) and kurtosis (peakedness of a dis- 
tribution) revealed significant differences, with the rmTBI 3d group 
having increased skewness (p = 0.026) and kurtosis (p = 0.034) com- 
pared to the Single group but only at the 14d post injury time point 
(Table 1). In contrast, the rmTBI 7d group had significant increases in 
skewness (p = 0.042) and kurtosis (p = 0.033) compared to Singles 
only at the 1 d post last injury time point (Table 1 ). 



3.2. Injury volume is dramatically increased at acute time points in 
rmTBI animals 

The anterior-posterior extent of the mTBI injury volume was assessed 
computationally and included all tissue abnormalities, including blood 
(hypo-intense) and edema (hyper-intense) (Fig. 3A). A Single mTBI 
resulted in an average total brain lesion volume of 1.4 ±0.002% of total 
brain volume (lesion volume approximately 15.7 ± 1.6 ul) at Id post 
first injury, which appeared to have a heterogeneous mixture of blood 
and edema (arrowheads; Fig. 3A). Following a second mTBI (Id post 
last) the rmTBI 3d group exhibited increased edematous tissue (white 
arrowheads) at the site of the second impact, while a more heteroge- 
neous lesion appeared at the site of the initial injuiy (Fig. 3A). However 
in rmTBI 7d animals, injury was composed primarily of blood at the site 
of the initial impact, while a larger injury of similar composition to that 
of a Single mTBI was observed at the site of second impact (white 
arrowheads; Fig. 3A). At the time of last MRI, little or no observable ab- 
normalities were detected within the tissues (Fig. 3A). Volumetric 3D re- 
construction of mTBI volumes illustrates the location and extent of the 
mild injury of our model revealing an increase in total lesion volume 
within rmTBI 7d animals compared to other groups (Fig. 3B). 



Temporal mTBI injury volumes were further assessed to deter- 
mine whether ongoing pathology from the first mTBI resulted in an 
exacerbation of lesion volume following a second mTBI (Fig. 3C). 
The rmTBI 7d group demonstrated a significantly increased lesion 
volume (p = 0.018) compared to other groups at Id post last injury 
(Fig. 3C). The observed increases in the rmTBI 7d lesion volumes 
were transient with lesion volumes similar to those observed in 
other groups by 14d post first injury (Fig. 3C). We then further 
dichotomized the mTBI injury volume by examining independently 
the first and second lesion volumes to test whether differences 
existed between the first and second injuries. This analysis revealed 
a significant increase (p = 0.029) in lesion volume at the site of the 
second injury within rmTBI 7d animals Id post last injury compared 
to Single animals, which was not observed in the rmTBI 3d group 
(Fig. 3D). The rmTBI injuries seen at the site of initial impact appeared 
similar to that observed in Single animals at Id post last injury 
(Fig. 3D). By 14d post injury (first: 14d post first, second: 14d post 
last injury), significantly decreased lesion volumes were observed in 
all rmTBI animals compared to Singles (Fig. 3D). 

3.3. Interval between rmTBI dictates injury characteristics: edema versus 
blood 

Computational analyses of the first and second lesions were per- 
formed to determine the relative contribution of blood and edema 
(based on T2 values) at each injury site (Fig. 4). At Id post last injury, 
computational analysis of the voxels containing blood (green), edema 
(red) and NAB (blue) at the injury sites clearly demonstrates that lesion 
within Single mTBI animals had a mixture of voxels containing blood, 
edema and NAB as was seen in the MRI (arrowheads; Fig. 4A). In contrast 
rmTBI 3d animals exhibited predominately edema (red) at the site of the 
second mTBI, as was observed in T2 MRI (hyper-intense signal; arrow- 
heads). The site of the first injury within rmTBI 7d animals exhibited al- 
most entirely blood (green) containing voxels, which can also be seen 
in T2 and SWI images (arrowheads; Fig. 4A). 

Quantification of these voxel groups following a second impact in 
rmTBI 3d animals demonstrated a significant (p = 0.038) increase in 
the number of edema voxels (increased volume) Id post last injury, 
which was dramatically resolved by 14d post injury (first: 14d post 
first, second: 14d post last injury) compared to Single animals 
(Fig. 4B). Further the rmTBI 3d group had a significantly (p = 0.003) 
greater edema volume than rmTBI 7d animals. The site of initial injury 
within the rmTBI 3d group also showed elevated edema volumes com- 
pared to other groups, though this was not significant (Fig. 4B). In con- 
trast, the second injury in rmTBI 7d animals exhibited similar edema 
volumes to that observed after a Single mTBI, while the first lesion 
had a significant (p = 0.001) decrease in edematous tissue compared 
to the Single group at Id post last injury (Fig. 4B). 

Both injury sites of rmTBI 7d animals demonstrated a significant in- 
crease (1st: p = 0.0001; 2nd: p = 0.025) in blood volume compared to 
a Single mTBI at Id post last injury (Fig. 4C). However, while there was 
an increase in the total amount of blood seen in the rmTBI 7d animals, 
the first and second blood volumes were not significantly different 
from each other. Further the rmTBI 7d group contained significantly 
more blood than the rmTBI 3d group (p = 0.04). The first lesion within 
rmTBI 3d animals had a similar blood volume to that seen after a Single 
mTBI Id post injuiy, while the site of second injuiy had a smaller blood 
volume compared to other injuiy groups (Fig. 4C). Significant differences 
(p = 0.047) in blood volume were still observed 14d following the initial 
injuiy, while blood deposition at the site of the second mTBI appeared 
similar between the rmTBI groups at 14d post second injury (Fig. 4C). 

Prussian blue staining for extravascular blood at the site of the lesions 
was performed 14d post first injury and supported our MRI analysis, 
demonstrating increased blood deposition within the rmTBI 7d animals 
compared to other groups (Fig. 4D). Blood deposition also appeared 
concentrated underneath the impact site within the Single and rmTBI 
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Fig. 3. rmTBI 7d animals have increased lesion volumes. A. Representative T2WIs illustrating the presence of abnormal tissue after the initial (Id post first; yellow arrows) and second (Id 
post last; white arrows) mTBl. B. 3D reconstruction of injury volumes illustrates the mild nature of the mTBI injury, where average total lesion volumes from Single (1.4%), rmTBI 3d (0.9%) 
and rmTBI 7d (2.0%) animals were collected Id post last injury. C.The temporal evolution of total (first + second) mTBl lesion volumes over the experimental period, demonstrates a sig- 
nificantly increased lesion volume in rmTBI 7d animals Id post last injury compared to other groups. D. Evaluation of the first and second injuries revealed a transient increase (p = 0.029) 
in the second lesion volume of rmTBI 7d animals compared to Singles Id post last injury. Normalized data presented as means ±SEM, where *p<0.05. 



3d groups, while the rmTBI 7d animals had a more diffuse deposition 
(see below, Fig. 4D). 

3.4. Activated microglia demonstrates spatial overlap with blood 
deposition 

We undertook an analysis of the distribution of Prussian blue 
staining (blood) at 14d post first injury (Fig. 5), where we observed a 
modest increase in the blood containing area at the site of the first 
impact in rmTBI 7d animals (Fig. 5A). In contrast, the site of the second 
injury in rmTBI 3d animals revealed blood deposition primarily under- 
neath the impact site, while rmTBI 7d animals demonstrated a dramatic 
increase in width and depth of blood (Fig. 5A), often down to the corpus 
callosum. 

Measurements of activated microglia (IBA1) containing regions 
revealed a similar profile with locations encompassing blood deposition 
(Fig. 5A). The site of initial injury demonstrated a smaller area of 
microglial activation compared to blood containing zones in both 
rmTBI 3d and Single groups, while regions were approximately the 
same size in rmTBI 7d animals (Fig. 5A). In contrast, the site of second 



injury in both rmTBI groups had an increased dispersion of microglial 
activation beyond the regions containing blood, which was most prom- 
inent within the rmTBI 7d group (Fig. 5A). Staining for astrocytes 
(GFAP) revealed increased reactivity only at the site of second injury 
within rmTBI groups, which was confined to the immediate vicinity of 
the impact site (data not shown). 

A modest number of Fluoro-Jade positive cell-bodies were observed 
at the site of the second impact in rmTBI 7d animals suggesting some 
ongoing neuronal degeneration, which was not seen in the first injury 
site or in other groups (data not shown). However in all mTBI groups 
punctate staining was observed, reminiscent of axonal or dendritic de- 
bris, at the site of the first injury but not in the left non-injured cortex 
of Single animals suggesting ongoing microstructural damage (data 
not shown). 

4. Discussion 

Using a novel computational approach we report the temporal de- 
tection and quantification of blood and edema for the first time, pro- 
viding new insight into the progression of neuropathology following 
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Fig. 4. Voxel-wise analysis reveals increased blood volume in the rmTBI 7d animals. A. MR1 taken Id post last injury reveals abnormal tissue on T2 (yellow arrows) and blood de- 
position on SWI (yellow arrows) following mTBI. Computationally color-coded images illustrate voxel characterization (green = blood, red = edema, blue = normal appearing 
brain) in these animals. B. Edema volume analysis within the first and second lesions revealed changes in edematous tissue at Id post last injury, while no differences were ob- 
served between groups at 14d post injury (first: 14d post first, second: 14d post last injury). C. Analysis of blood volume within the first and second lesions demonstrated signif- 
icantly increased blood deposition within rmTBI 7d animals at Id post last injury and between the first lesions at 14d post injury (first: 14d post first, second: 14d post last injury). 
D. Prussian blue staining in the first and second lesions of animals 14d post first injury at the site of maximal lesion shows increased blood deposition in the tissues from rmTBI 7d 
animals. Normalized graphs presented as means ±SEM where *p<0.05; cal bar=100 um. 



a unique model of rmTBI. Our novel findings include: 1 ) T2 histogram 
evidence of increased edematous tissue in rmTBI 3d animals and in- 
creased blood deposition in rmTBI 7d animals. 2) Increased lesion vol- 
ume within rmTBI 7d animals, which was not observed in the rmTBI 
3d animals. 3) Hemorrhagic progression of the lesion was detected, 
using a voxel-wise comparison, within the lesion of rmTBI 7d animals, 
while increased edematous tissue was observed in rmTBI 3d animals. 

4) Prussian blue staining confirmed the presence and revealed the 
spatial distribution of blood within the cortex 14d post first injury. 

5) At the site of the second injury, an increased inflammatory re- 
sponse was detected in both rmTBI groups, while ongoing neuronal 
death was only observed in rmTBI 7d animals. Taken together, our 
study further demonstrates that there are differential increases in tis- 
sue vulnerability following rmTBI and that the underlying pathology 
is dependent upon the interval between and location of mTBI events. 

4.1. Lesion histograms predict temporal changes in rmTBI pathology 

Initially, we used voxel-wise T2 histograms to analyze rmTBI le- 
sions to evaluate global temporal changes in tissue pathology. This 
is the first report using T2-histograms in an experimental model of 
mTBI, as there are no previous publications to our knowledge. In con- 
trast, other clinical studies and animal models of central nervous sys- 
tem (CNS) injuries including tumors, spinal cord injury and ischemia, 



have used histogram analysis (Emblem et al., 2008; Kopelman et al., 
2005; Nevo et al., 2001; Veltkamp et al., 2005). 

Clinically, histograms have also been employed to quickly assess 
whole brain and tissue-specific (white and gray matter) MRI changes 
in mild to severe TBI (Lipton et al„ 2008; Marquez de la Plata et al., 
2011; Shanmuganathan et al„ 2004). While histogram analysis has 
demonstrated some ambiguity in sensitivity from clinical data others 
have reported its effectiveness to discriminate TBI. Whole-brain histo- 
grams in mTBI patients undergoing diffusion tensor imaging (DTI) 
found no differences in the mean diffusivity (MD) and fractional anisot- 
ropy (FA) measurements between injured and control patients (Inglese 
et al., 2005). However when the histogram analysis was focused to ROIs 
drawn within white matter tracts susceptible to mTBI, increased MD 
and decreased FA were observed in injured patients (Inglese et al., 
2005). In contrast, application of global histograms (whole white mat- 
ter) from DTI images were able to identify mTBI, reporting decreased 
FA without the need for ROIs (Benson et al., 2007). These clinical studies 
suggest that histograms have the potential to detect subtle tissue 
changes associated with mild injuries when analysis is focused on spe- 
cific brain regions rather than performing whole-brain histogram anal- 
ysis. In our study we observed temporal histogram skewness and 
kurtosis changes within defined mTBI lesions, suggesting tissue hetero- 
geneity associated with edema resolution and ongoing blood deposi- 
tion. These findings confirm that histogram analysis has the potential 



Fig. 5. Microglial activation and blood deposition demonstrate spatial overlap within the lesion sites. A. A depiction of the rat cortex illustrates the locations of observable tissue damage 
where blood and microglial measurements were taken. Depth and width graphs demonstrate a spatial overlap in the blood and activated microglial localization at 14d post first injury. 
Single data in the second injury (dotted line) is from the injured cortex (first injury) for comparison. B. Representative Prussian blue and IBA1 images, taken from adjacent sections of 
rmTBI 7d (second injury) and Single (first injury) animals, show the localization of blood and activated microglia within the cortices, cal bar = 100 um; inset cal bar = 20 um. 
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to discriminate between the mTBI groups as the lesions progress and re- 
solve over time. Our results set a foundation for future comparisons of 
whole brain and regional histogram analysis so that additional methods 
can be developed for efficient and accurate comparisons of clinically de- 
fined mTBI. 

4.2. Exacerbation of rmlBl lesion volume is dependent on the time 
interval between injuries 

We found that tissue damage was exacerbated following a second 
mTBI, particularly, when the injuries were 7d apart. The increased le- 
sion volume was predominantly due to exacerbated tissue damage at 
the site of the second injury. This suggests that a Single mTBI can af- 
fect distant brain regions in such a way that they remain vulnerable to 
a subsequent injury. 

The neuropathology and long-term consequences of rmTBl are not 
well understood. Rodent models of experimental rmTBI have demon- 
strated that there exists a window of vulnerability in which two inju- 
ries can result in exacerbated tissue damage (Longhi et al., 2005; 
Vagnozzi et al., 2007). Animals sustaining mTBls 3 or 5 days apart 
on the same hemisphere exhibited increased cognitive deficits for 
up to a week following injury and cell death and axonal disruption 
72 h following the last injury (Longhi et al., 2005). Similarly, an ex- 
amination of changes in energy metabolism using MR spectroscopy, 
demonstrated the greatest changes following mTBIs sustained 3d 
apart and which persisted for up to a week (Vagnozzi et al., 2007). 
However, the changes observed by Longhi et al. (2005) were not ob- 
served in animals receiving two impacts 7 days apart, while alter- 
ations in metabolism described by Vagnozzi et al. (2007) were not 
observed in animals sustaining injuries 5 days apart. Collectively, 
these studies suggest that by > 5 days the tissue has had enough 
time to heal and so a second impact was similar to sustaining a Single 
mTBI when delivered to the same hemisphere. Additional rmTBI 
studies inducing two impacts 24 h apart have further supported 
these findings, where in the absence of gross histological abnormali- 
ties, increased neurological deficits were persistent for several 
weeks with an acute increase in axonal injury and increased inflam- 
mation (Laurer et al., 2001; Shitaka et al., 2011). 

However these models induced rmTBI to the same anatomical lo- 
cation, while our model induces mTBI injuries distant from each 
other. This suggests that brain vulnerability is dependent on not 
only the interval between injuries but also the location where the in- 
juries are sustained. Patients sustaining an mTBI demonstrated meta- 
bolic tissue deficits for approximately 1 month following injury while 
those sustaining a second mTBI within 10-13 days of the first demon- 
strated a prolonged alteration (Vagnozzi et al., 2008). While perhaps 
brain regions close to each other may have a shorter window of vul- 
nerability, locations farther apart appear to have a delayed or length- 
ened period of vulnerability. 

4.3. HPC primarily contributes to the lesion volume of rmTBI 7d animals 

Computational analysis of blood and edema contribution to the ab- 
normal tissue volume revealed that rmTBI 3d animals had increased 
edema volumes, while the rmTBI 7d group exhibited increased blood 
deposition. By 14d post injury (14d post first and 14d post last injury) 
the lesions of both rmTBI groups appeared to have a persistent increase 
in blood deposition, while edema resolved. It is important to note that 
extravascular blood can demonstrate a temporal change in MRI signa- 
ture, as observed in an animal model of intracranial hemorrhage (ICH) 
(Belayev et al., 2007). However due to the mild nature of our TBI 
model and the ongoing blood deposition, as a result of hemorrhagic 
progression, we would expect a different temporal course of extravas- 
cular blood deposition than that observed following ICH. 

Our findings of increased blood deposition with time concurs with 
the concept of hemorrhagic progression of the contusion (HPC), 



which has been modeled experimentally in moderate to severe trau- 
matic brain and spinal cord injury (Gerzanich et al., 2009; Simard 
et al., 2009). Within 12 h of a moderate to severe TBI, HPC was observed 
where the blood volume expanded (2 fold) into the hippocampus and 
thalamus (Simard et al., 2009). Similarly, in our rmTBI 7d model, at Id 
post last injury we observed a 2.5 fold increase in blood volume at the 
site of the initial injury in the rmTBI 7d group. Several studies have dem- 
onstrated concordance between Prussian blue staining of tissue iron de- 
position and T2 weighted MRI (Belayev et al., 2007; Gupta et al., 2012; 
Thulborn et al., 1990). Our histological staining for extravascular blood 
further supported the expansion of bleeding in the rmTBI 7d animals, 
which was not observed in other mTBI groups nor has been reported 
in other models of rmTBI. 

Models of CNS injury have identified the sulfonylurea receptorl- 
regulated NC Ca _ ATP channel (SUR1/TRPM4) as a key mediator in the 
development of HPC (Gerzanich et al., 2009; Simard et al., 2009). 
These molecules are upregulated following trauma and upon deple- 
tion of ATP, a deregulation of this channel occurs leading to capillary 
fragmentation and hemorrhage (Simard et al., 2010). HPC in our 
model was identified using temporal MRI at acute time points and 
Prussian blue staining at sub-acute times following injury. Our tissue 
was collected at late time points, thus we could not identify the cellu- 
lar cascades, which occur within hours to days following trauma 
(Simard et al., 2009). 

Clinically, HPC has been identified using CT, where approximately 
45-50% of patients develop this pathology following a mild to severe 
TBI (Alahmadi et al., 2010; Narayan et al., 2008; Sifri et al., 2004, 
2006). A clinical study quantifying changes in hematoma size, using 
CT following a TBI (mild to severe), found that 4% shrank, 58% were 
unchanged and 38% increased in volume within 72 h (Chang et al., 
2006). Clinical evidence of HPC suggests that medical imaging can be 
used to identify this pathology, though future experimental studies of 
the underlying cascades need to be performed to validate and better un- 
derstand the underlying HPC mechanisms. 

Another potential mechanism for our observed increase in blood 
volume could be due to extracellular matrix remodeling and subse- 
quent revascularization of the tissue. In models of ischemic stroke, 
hemorrhagic transformation has been observed and associated with 
breakdown of the extracellular matrix (ECM), which contributes to mi- 
crovascular integrity (Fukuda et al., 2004; Hamann et al., 1995). 
Metalloproteinases (MMP) have been shown to be associated with 
microvessel ECM degradation and blood brain barrier permeability fol- 
lowing experimental models of stroke (Fukuda et al., 2004; Gidday et 
al., 2005). In models of moderate to severe TBI, MMPs are upregulated 
with maximal levels being reached at 24 h and can remain high for up 
to 7 days following injury (Hayashi et al., 2009; Wang et al., 2000). In 
cell culture models of ischemia, microglia increased secretion of the in- 
active form of MMP-9, which in non-human primate models of stroke 
demonstrated a linear relationship with blood deposition (Del Zoppo 
et al., 2012). Treatment of moderate TBI with Minocycline, a drug 
shown to suppress microglial activation, inhibited MMPs and resulted 
in increased synaptic preservation (Ding et al., 2009). Additionally, 
knockout of MMP-9 demonstrated decreased lesion volumes and im- 
proved neurological outcomes following moderate TBI compared to 
control mice (Wang et al., 2000). Future examination could include tis- 
sue collected at these early time points for changes in the ECM as a re- 
sult of MMPs. Collectively these studies demonstrate that MMPs play 
a role in tissue exacerbation following trauma though their role in 
HPC is not clear. 

4.4. Increased neuroinflammation at the site of the second injury in 
rmTBI animals 

Extravascular blood within brain tissues can result in apoptosis and 
initiate a neuroinflammatory response, including peripheral leukocyte 
infiltration and microglial activation (Gong et al., 2000; Xue and Del 
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Bigio, 2000). We found that there was a distinct overlap of brain regions 
containing activated microglia and extravascular blood, although 
microglial responses encompassed a slightly larger area at the site of 
the second impact, particularly in the rmTBl 7d group. Microglia have 
been identified as potential aggravators of CNS tissue damage. During 
CNS infection or injury classically activated microglia release reactive 
oxygen species causing further tissue damage, but suppression of mi- 
croglia demonstrates decreased TBI related tissue loss (Boje and Arora, 
1992; Homsi et al., 2009). A study of the inflammatory response to 
microhemorrhages within brain tissue found that microglia within 
200 urn of blood can remain responsive for up to 7 days (Rosidi et al., 
201 1 ). Our observed increase in microglial response and its detrimental 
effects is consistent with increased lesion volumes observed at the site 
of second injury within rmTBl 7d animals. 

4.5. Conclusions 

Our new model of rmTBl, in which injuries were induced to opposite 
hemispheres, demonstrated that there exists a window of tissue vulner- 
ability following an initial mTBI. We previously hypothesized that ani- 
mals receiving impacts 7d apart would not exhibit increased tissue 
damage, while injuries sustained 3d apart would result in larger lesion 
volumes specifically at the site of the second injuiy. However, our data 
supports that this temporal window occurs at a later time point than pre- 
viously described (7d vs. 3d) and is likely the result of the second injuiy 
location being distant from the first. Analysis of cortical abnormalities 
using our rapid computational approach revealed increased blood depo- 
sition within rmTBl 7d animals consistent with HPC, which was not ob- 
served in the rmTBl 3d group. This increased blood deposition was also 
observed to co-localize with regions of activated microglia, which may 
further cause tissue damage at regions distant from the impact site lead- 
ing to the increased lesion volumes observed in rmTBl 7d animals. The 
vulnerability to a second impact resulted in both a greater volume of tis- 
sue damage and increased blood deposition following a second tempo- 
rally and spatially distant mTBI. The differences in observed blood and 
edema volumes between the rmTBl groups are most likely the result of 
temporal cellular cascades, such as HPC, initiated by an initial injuiy. 
However additional investigations are needed to uncover the mecha- 
nisms leading to increased brain vulnerability and blood deposition 
within rmTBl 7d animals and determine the long-term effects of this ob- 
served tissue exacerbation. Further computational analysis methods 
need to be developed for experimental models of TBI so that they can 
be optimized for future clinical use. 
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